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Abstract. We study the problem of continuum percolation in infinite volume 
Gibbs measures for particles with an attractive pair potential, with a focus 
on low temperatures (large The main results are bounds on percolation 
thresholds p±(/3) in terms of the density rather than the chemical potential 
or activity. In addition, we prove a variational formula for a large deviations 
rate function for cluster size distributions. This formula establishes a link with 
the Gibbs variational principle and a form of equivalence of ensembles, and 
allows us to combine knowledge on finite volume, canonical Gibbs measures 
with infinite volume, grand-canonical Gibbs measures. 

Keywords: continuum percolation, stochastic geometry, point processes, large 
deviations; Gibbs measures, Gibbs variational principle. 

MSG subject classification: 60D05, 60K35, 82C43. 



1. Introduction 

The present article is concerned with percolation properties for Gibbsian point 
processes. We arc interested in infinite volume Gibbs measures (in the sense of the 
Dobrushin-Lanford-Rucllc conditions) for particles in M. d interacting via an attrac- 
tive, finite range pair potential. The dimension is two or higher, d > 2. Around 
each particle a; of a random configuration cj, draw a ball of radius R, for some fixed 
R > 0. Percolation occurs if the region in M. d covered by the union of such balls, 
U xe . u B(x,R), has an unbounded connected component, with positive probability. 
In the notation of Meester and Roy [MR96] , our problem is a Boolean percolation 
model (X, p) driven by a Gibbsian point process X and with deterministic radius 
p = R. 

This problem has been studied before [M75l fZ08| fP Y09) lArT2] . for both repulsive 
pair potentials and potentials with an attractive part. Miirmann |M75| investigated 
finite-range potentials and gave a sufficient condition for the absence of percolation; 
a different proof, with an extension to tempered boundary conditions for attractive 
potentials, was given by Zessin |Z08j . His proof builds on integration by parts 
for Gibbsian measures. Pechersky and Yambartscv [PY09J proved a criterion for 
absence of percolation with the help of a coupled branching process; their result 
does not require that the potential has compact support. In addition, Pechersky 
and Yambartsev gave a sufficient condition for the presence of percolation, valid 
in dimension 2, for attractive potentials with possibly unbounded support. An 
analogous result, for hard spheres in dimension 2 was shown by Aristoff [Arl2] , 

The cited works all formulate criteria in terms of the activity z - i.e., the intensity 
parameter of some a priori Poisson point process - rather than the density p, 
which for interacting particles is a non-trivial function of z. As a consequence, for 
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attractive potentials, the cited results cannot distinguish between two very different 
physical pictures. First, percolation might be a high-density effect, as expected for 
hard spheres; second, it might be an energetic effect - at low temperature, it may 
happen that the interaction favors the formation of large connected components, 
which because of entropy may coexist with large almost empty regions of space 
containing a few small components. In such a situation the density threshold for 
percolation may be very small, as suggested by recent results jJKMllj on large 
deviations for cluster size distributions in the canonical ensemble. 

The aim of the present article is, therefore, to give bounds on percolation and 
non-percolation thresholds in terms of the density rather than the activity. Our 
main result is Theorem 13.61 below, valid for finite-range, attractive potentials. 
Quickly summarized, Theorem 13.61 states that there are curves p±(/3) such that 
if Pp lP is a shift-invariant Gibbs measure at density p and inverse temperature j3, 
the following holds. 

• If p < p~{(3), there are only bounded connected components, P^ p -almost 
surely; we have P-(f3) = exp(— (3v*(l + o(l))) for suitable v* > 0. 

• If p > p + (/3), there is an unbounded connected component, Pr p - almost 
surely; as f3 — > oo, p+((3) —> po > for suitable po > 0. 

We expect that for very large (3 and P-{(3) < p < p+(ft), there are non-ergodic 
Gibbs measures with percolation probability strictly between and 1. We have 
no proof of this conjecture; sec, however, Proposition 13.71 for some preliminary 
evidence, and Appendix [S] for relevant lattice gas results. 

The key technical tool for the proof of Proposition 13.71 is a variational formula 
for a large deviations rate for cluster size distributions in the canonical ensemble 
(Theorem 13. 2|) : the large deviations principle was investigated in JKM11 . The 
variational formula establishes a relation with the Gibbs variational principle and 
allows us to apply a form of equivalence of ensemble |G95] . This in turn enables us 
to combine the knowledge for grand-canonical, infinite volume Gibbs measures in 
[M751IPY09] with results on the canonical ensemble jJKMllj . 

We conclude the introduction with a word of caution on the physical interpre- 
tation of percolation. We should stress that a percolation transition need not be a 
phase transition - for the ideal gas at activity z (Poisson point process with intensity 
z), the pressure is an analytic function no matter the value of z, even though there 
is a percolation transition at high enough z |MR96j . Nevertheless, for attractive 
pair potentials and at low temperature, the percolation transition might coincide 
with a phase transition. In fact, for nearest neighbor attractive lattice gases (or 
Ising model), and temperatures below some threshold T + , the percolation transi- 
tion and the phase transition coincide (see the review by Georgii, Haggstrom and 
Maes [GHMOlj ). In dimension two, T + equals the Curie temperature Tc, but in 
higher dimensions T + < Tc [ABL87] , illustrating again that percolation and phase 
transition in general do not coincide. For the reader's convenience, we summarize 
some relevant results in Appendix [A] 

The remainder of the article is organized as follows. In Sections [5] and [3] we 
formulate the setting and our results. Section [4] defines the topology of local con- 
vergence and summarizes continuity properties of important functions such as the 
relative entropy rate. Sections [5] to [7] arc devoted to the proofs. 
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2. Setting 

2.1. Pair potential. The pair potential is a function v : [0,oo) -> EU {oo} that 
serves to define the total energy of an A-particle configuration 

U N (xi,...,x N ) := ^ v {\ x i- x j\)-> x l-> • ■ • ' x n G Rd - 

l<i<j<N 

By a slight abuse of notation, we shall drop the subscript and write U (x) instead 
of Un{x). 

Assumption 1. The pair potential satisfies the following basic assumptions: 

• Either v is everywhere finite or there is a r^ > such that v(r) = oo for 
r < rhc and v(r) < oo for r > rh c - (We impose no condition on v(r\ lc ).) 

• v has compact support: n := supsuppu < oo. 

• v is bounded from below: inf v > — oo. 

• v has an attractive tail: for suitable ro < ri and all r G (ro, n), u(r) < 0. 

If rjj C > 0, we say that v has a /iard core. If u has compact support, we shall also 
say that v has finite range. For k € Z d , let C(k) be the unit cube [ki + 1) x • • • x 
[fc rf , + 1) and 

7V c(fc) (xi,...,z N ) := \{j G {l,...,JV}|a?j £ C(fe)}| 
the number of particles in C(k). 

Assumption 2. The pair potential v is super stable: there are constants a > 0, 
& < oo such that for all N G N and all x = (xi, . . . ,xjv) G (K d ) Ar , 

> ^ (a7V c(fc) (x) 2 - bN c{k) (xfj . (1) 
feez d 

Every supcrstablc interaction is stable, since Eq. ([1]) implies in particular that 
U(x\, . . . , xn) > —bN. If u satisfies Assumption [T] and is non-intcgrably divergent 
at the origin, then f is supcrstable; see |R70| IR69j for a proof and other sufficient 
conditions. 

Assumption 3. The potential is integrable in {v < oo}: Jj x i >rh ^(l^l) |dx < oo. 

Assumption 4. There is an r m ; n > such that for every N efj, U has a minimizcr 
(xi, . . . , a; at) G (IR c ') Ar with interparticle distance bounded from below by r m ; n : 

1 < i < j < N => \xi - Xj | > r min . 

In addition, v is Holder-continuous in [r m i n ,oo). 

A sufficient condition for the lower bound on interparticle distances is that 
v(r)/r d — > oo as r —> 0, as can be shown along }T06| Lemma 2.2]. 

Assumption 5. There is a C > such that for every G N, U has a minimizcr 
{xx, . . . ,x N ) G (M d ) Ar with diameter bounded by CA 1 ^: 

maxjla;, - Xj\ | 1 < i < j < N} < CN 1/d . 

Under Assumption [T] this condition is trivially fulfilled in dimension 1. In di- 
mension 2, sufficient conditions are given, for example, in |T06j . where much more 
is proven on the ground states. To the best of our knowledge, there is no result in 
dimension 3 or higher; in fact, providing upper bounds on interparticle distances 
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for Lennard- Jones type interactions seems to be a non-trivial problem in non-linear 
optimization, see the article by Blanc |B04| and the references therein. 

Let us briefly comment on our conditions on the pair potential. Assumption [T] 
simply defines the class of pair potentials we are interested in. Supcrstability as in 
Assumption [2] is a standard condition that ensures the existence of infinite volume 
Gibbs measures, see the next subsection. The intcgrability assumption |3] will allow 
us to use a bound on Mayer expansions going back to Brydes and Fcderbush |BF78| . 
Assumptions 0] and [5] are needed for precise statements about Gibbs measures at 
low temperature /3 _1 and densities above some threshold of the form exp(— f3v*). 

2.2. Infinite volume Gibbs measures. Let £1 be the set of locally finite point 
configurations, 

ft := {u C R d | Vr > : \uj fl B(0, r)\ < oo}, 

with P(0,r) the open ball of radius r centered at the origin. We equip fi with 
the er-algebra T generated by the counting variables Nb(uj) := \oj fl B\, B C M d 
Borel-measurable, and denote the probability measures on (f2, IF) with the letter 
V . The following subsets of V will be relevant for us: shift-invariant measures Vg, 
tempered measures, and infinite- volume Gibbs measures £/(/3, fi); we proceed with 
their definition. 

For x G K d , the shift 6 X : Q -> Q is defined by 9 x u := {y — x \ y G oj}. A 
measure P G V is shift-invariant if P(8 X (A)) = P(A) for all measurable A C ft and 
all x G R d . The collection of shift-invariant measures is denoted Vg. We say that 
P G V is tempered if for P-almost all w, 

3t(w) > W G N : ^ (^c(fc)H) * < (2) 

kez d n[-e,e] d 

where C(k) is the unit cube [ki, ki + 1) x ■ • • x fed + 1). Fix (3 > and /i G M. 
Let A C M d be a bounded Borel set and ( g O a configuration satisfying the 
tempcredness condition ([2]). For n G N, define the measure Q„ on A" as the 
measure with Lebesgue density 

f n 
— exp(-/3 ^ u(|a* -Xjl) v (\ x ^y\))- 

l<i<j<n i=l yGCnA c 

The image of Q n under (xi, . . . , x n ) <-> {xi, . . . , x n } is a measure on O which we 
denote again by Q n . Let Qo be the probability measure on fl which gives probability 
1 to the event that u = 0. Define p£ l(ttl A|c by 

^ oo 

The normalization Sa|£ (/?,£«) is defined by the requirement that P/3,^ a|£ is a prob- 
ability measure on O. 

We say that P G P is a (/3, n) -Gibbs measure if and only if P is tempered and 
for every A and every measurable / : Q — > [0, oo), 

\ P(dw)/(w)= / P(dC) / P^, M , A | C (d W )/(u;U(CnA c )). (3) 

We denote the set of (/3, /i)-Gibbs measures as £(/3, /x); the sets Q(fl,(j.), G 
M+xK, are pairwise disjoint [G95| Remark 3.7]. The existence of Gibbs measures for 
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lower regular, superstable potentials - this class includes our potentials satisfying 
Assumptions Q] and [2] - is shown in [R70| lD70) : for all (3 > 0, fj, e R, 
and Q(f3,p.) C\Ve are non-empty, convex sets. With suitable topologies, they are 
Choquet simpliccs and, in particular, compact |R70[ Theorems 5.6 and 5.8]. 

2.3. Palm measure; energy and entropy densities. For each P e P«, there 
is a unique finite measure P° on (f2, !F) such that for all measurable, non-negative 
functions / : R d x fl -> [0, oo), 

/ Y j .f{x,e x uj)P{Auj)= [ f f(x,u)dxP°(du). (4) 
Jn X&UJ Jn Jr* 

P° is the Palm measure of P |DVJ081 Chapter 13]. It is known that P°(0 £ w) = 
and P°(Q) = p(P), where p(P) is the expected number of particles in [0, l) d , called 
density. The defining equation (j4j, specialized to f(x, w) = lc(x)g(ui) with C C K d 
Borel-measurable, yields the useful identity 

/ <?MP°(d w ) = jL f g{e x uj)P{Au). (5) 

|C| is the Lebesgue volume of C. 

Remark. The Palm measure is the continuum analogue of a simple lattice object: 
if O = {uj | u> C Z d } = {0, 1} Z and P is a shift-invariant measure on f2, the Palm 
measure becomes P°(A) = P(A n {0 e w}). 

For a shift-invariant, tempered measure P, we define the expected energy per 
unit volume, or energy density, as 

u ( p ) := l I (E y (iyi)) po ( dw )^ / ( E $>(|y-*i))p(cM. 

Ji yew ^ n x6wn[o,i] d yew 

i7(P) takes values in RU {oo}. For superstable potentials as in Eq. ([!}, the energy 
is bounded from below: U(P) > —b 2 /Aa. 

The entropy per unit volume, or entropy density is defined as a relative entropy 
rate. Let Q G Ve be the Poisson point process with intensity 1 and P € Vg. For 
A = [— L, L] d C M d , let Pa be the image of P under the projection map wnjuflA. 
Define Qa in a similar way. The relative entropy is 

I(Pa;Qa) ■= I /(w)log/(w)Q A (dw)= [ log f(cj)P A (dw) 
Jn Jn 

if Pa has Radon-Nikodym derivative dP\/dQA = /, and I(P\;Qa) ■= oo if Pa s 
not absolutely continuous with respect to Q\. The limit 

5(P):=1- lim i-/(P A ;Q A ) eKU{-oo} (6) 

|A|->oo A| 

exists for all P £ Ve, see |GZ93| . Continuity properties of U(P) and S(P) with 
respect to a suitable topology on V are recalled in Section 2] below. 

Remark. The additive constant 1 is included in Eq. [S] for aesthetic reasons; if we 
did not include it, we would need an additive constant in Eq. (|9]). 
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. e -«|A| 



P*,a(w) := ^ l(|C wn A(i)| = k )- 



Example. Let P be a Poisson point process with intensity parameter z. Then 

QA) = £ i l0g |A|«exp(-|A|) J 
= (l-z)|A| + z|A|logz 
thus 5(P) = -z(logz - 1). 

2.4. Cluster densities. Let ri > be the range of the potential as in Assumption 
[T] Fix R > n. For w S fi, let G w be the graph with vertex set ui and edge set 
{{x, y} | x,y e w, < |x — y| < i?}. When a; € w, let C u (x) C w be the connected 
component of x in G u . When x £ M. d \u>, we set C UJ (x) := and |C w (x)| := 0. For 
P (z Vb and fcgff, the expected number of k-clusters per unit volume or k-cluster 
density is 

Pfc (P):=fc- 1 P°(|C(0)|=fc) = i- / £ l(\C u (x)\=k)P(duj). (7) 

The last identity holds for every Borel set C, compare Eq. <j5j) . Note that for 
every P £ Vg, X^fcli kpk{P) < p(-P)- For later purpose we also define finite volume 
empirical densities as 

f 

fclAI 

1 1 iGu 

Thus /9fc,A(w) is the number of fc-clusters in wflA, divided by the volume |A|. We 
have, for all lu G O, Sfi^ kp kA (uj) = AT A (w)/|A|. 

2.5. Large deviations for cluster size distributions. Finally we recall a large 
deviation principle shown in jJKMllj . Fix R > n. For A = [0, L] d and N e N, let 

Ziifi, N) := / e-W^'-'^dxi ■■■dx N . 

iV! J A N 

be the canonical partition function and P^ jv A the probability measure on A N with 
density N\~ 1 Za((3, N)^ 1 exp(—(3U(x 1 , . . . ,xn)). The image of P/3,at,a under the 
mapping R w — > 17, (xi, . . . , xjv) M- {xi, . . . , x^r} is a measure on (Q, J 7 ), for which 
by a slight abuse of notation we use the same letter P^ Aj at. 

Equip with the product topology and the associated Borel a-algebra, and let 
p A : (CI, J 7 , P^jv.a) — > R+ be the random variable Pa(w) := (/9fc, A (w))fceN- We are 
interested in the behavior of p A in the thermodynamic limit 

N 

N -> oo, L -> oo, — -> p (8) 
|A| 

for p > 0. First we recall the definition of the free energy f(/3,p) and the close- 
packing density p cp . Consider the limit 

/(p» :=-lm-LlogZ A (p\JV) 

along ([8]). It is well-known |R69| that the limit exists and is finite when the density 
p is strictly smaller than the close-packing density p cp > 0, and is infinite when 
p > p cp . When the potential has no hard core (r^c = 0), we have p cp = oo and the 
limit f(/3,p) is finite for all p > 0. 
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The following holds jJKMllj : for all /? > and p g (0,p cp ), in the limit ©, 
the random variable Pa( w ) satisfies a large deviation principle with speed /3|A|. 
The rate function is of the form p, (pk)ken) — f{P, p) for a function f(f3, p, ■) : 
[0,oo) N -> RU{oo} that is convex and lower semi-continuous with compact sublevel 
sets. Moreover, 



Note that the rate function f(/3,p, •) (but not the free energy f((3,p)\) depends 
on the connectivity radius R > n; to lighten notation, however, we leave the R- 
dcpendcnce implicit. 



Here we formulate our main results. Section 13.11 provides a variational char- 
acterization of percolation. Sections 13.21 and 13.31 formulate bounds on percolation 
thresholds in terms of the chemical potential p = fi~ l logz and the density p. Sec- 
tion 13.11 on the one hand and Sections 13.21 and 13.31 on the other hand are logically 
independent, except for Proposition 13.71 

Throughout the remainder of the article we shall assume without further mention 
that the pair potential satisfies Assumptions [T] and [2j and the connectivity radius 
R is larger or equal to the potential range r\ . 

3.1. Variational characterization of percolation. Our first result expresses 
the rate function f((3,p, {pk)) of |JKM11] in terms of a variational problem. 

Theorem 3.1. For every p £ (0,p cp ), (3 > and (pk)k£fi S R+, 



f(P,P, {Pk)ken) 

= min{U(P) - /T^P) \PeV e , p{P) = p, Vfc e N : p k (P) = p k ) (9) 



with the convention min0 = oo. 

We recall Theorem 3.4 from [U95] . Fix j3 > and p g R. Then P g V e is 
a (j3, /u)-Gibbs measure if and only if it minimizes U(P) — /3^ 1 S(P) — pp(P); the 
minimum —p(/3, p) is minus the pressure. This is the Gibbs variational principle. 
Together with Theorem 13. 1[ it allows us to establish the following relation between 
minimizers of f(/3, p, •) and infinite volume Gibbs measures^ 

Theorem 3.2. Fix p g (0, p cp ) and (3 > 0. Then, for every (pfe)fcgN € R+, the 
following two statements are equivalent: 

(1) (p k )keN is a minimizer of f(fi,p,-)- 

(2) There is a chemical potential p € R and a shift-invariant Gibbs measure 
P eVed G(fi,p) such that p(P) = p and for all k g N, Pfe(-P) = /Ofe. 

The next elementary proposition establishes a relation between cluster densities 
and percolation, valid for general shift-invariant point processes. 

Proposition 3.3. Let P g Ve- The following statements are equivalent: 

^Georgii states his result under the additional assumption that the potential is non-intcgrably 
divergent at the origin. A close inspection of the proof shows, however, that we can dispense with 
this condition because our potentials have finite range; see the comment in IG95I before Lemma 



oo 




fc=l 



3. Results 
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(1) EZik P k(P)<p(P). 

(2) P{there is a cluster with infinitely many particles) > 0. 

(3) P {there is a cluster with infinite diameter) > 0. 

A quantitative relation is given in Eq. ([20]) below. We shall refer to both (2) and 
(3) as P(there is an infinite cluster) > 0, or more briefly as percolation. An imme- 
diate consequence of Theorem 13.21 and Prop. 13.31 is the following characterization 
of percolation in shift-invariant Gibbs measures. Set 

9e{fi, P ) := {Pe|J g(fi, n) | P G V e , P {P) = p}. 

jueK 

The Gibbs variational principle implies that Qe{f3,p) consists of the minimizers of 
U{P) — S{P) under the constraint p{P) = p; compare with the proof of Theorem 
EH below. 

Corollary 3.4. Let (3 > and p G (0, /o cp ) . Consider the following statements: 

(1) There is a P G Qg{fi,p) such that P{there is an infinite cluster) > 0. 

(2) f{P,p, ■) has a minimizer {pk) such that J2kLi ^Pk = P- 

(3) For every P G Qg{j3,p), P{there is an infinite cluster) = 1. 

(4) Every minimizer {pk) of f{/3, p, ■) satisfies J2kLi ^Pk < P- 
We have (1)^(2) and (3)^(4). 

A similar characterization holds for non-percolation. We note that Corollary 
13.41 establishes a relation between percolation for infinite volume Gibbs measures 
and cluster size distributions in finite volume canonical ensembles as examined in 
[JKMllj . 

3.2. Percolation thresholds: grand- canonical ensemble. Set E\ := and 

E 

E k := inf U{x 1 ,...,x k ) {k G N), := inf — - > -oo. 

(K^)fc ken k 

Let ro < r\ as in Assumption [T] and suppose that v is continuous in (ro,ri). Set 
—M := inf r>ro v{r) < 0. For m < M, choose r m G {ro,r\) such that v{r m ) < — m, 
see Figure 2 in jP Y09] . Note eoo < —dm. Consider the conditions 

V> > p+ VP G g{/3, p) P(there is an infinite P-cluster) = 1 (10) 

Vp < p- VP G Gifi,n) P(therc is an infinite P-cluster) = 0. (11) 

Set 

p+{/3-R) := inf{/i+ G K | p+ satisfies (fTOj) }. 
/i_(/3;P) := sup{/i_ gt|/i_ satisfies flT])}. 
Clearly R) < p + {/3; R). 

Theorem 3.5. Let R > r x , m G (0,M) and R m > \/d + 3 f m > Vrf + 3r . 
(1) Suppose that v satisfies Assumption^ Then 

e<x> < liminf p~{/3] R). 

In addition, for every p < and sufficiently large f3, there is a unique 
(/?, p)-Gibbs measure P; it is shift- invariant, has no infinite cluster (P- 
almost surely), and satisfies 

P°(|C,(0)| = k) = k Pk {P) < fce fe |P(0,P)| fc - 1 exp(-/3A : (e oo - p)). (12) 
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(2) Suppose that v is continuous in (ro,ri). Then 

limSUp /i+(/3; Rrn) < — 777 . 
fj— >oo 

We conjecture that for every fixed R > r%, under suitable conditions on v, we 
have P-(/3; i?) = p+((3; R) for sufficiently large /3 and 

lim yU_(/3;i?) = lim p + (/3;R) = e x . 

/3— >-oo /3— s-oo 

See Appendix lAl for the corresponding lattice gas result. 

3.3. Percolation thresholds: canonical ensemble. Let v* := inffc e pj(-£'fe — 

keoo). It is known that for potentials with an attractive tail, v* > |JKM11) . 
Consider the conditions 

Vp > p+yP € 0<?(/3, p) P(there is an infinite i?-cluster) = 1, (13) 

Vp < VP G £ e (/3, p) P(there is an infinite i?-clustcr) = 0. (14) 

Set 

p + (j3\R) := inf{p + e (0, p cp ) \ p + satisfies (HI])}, 

P-(/3]R) := sup{p_ € (0,/9 cp ) | P- satisfies (fTIj) }. 

Clearly p_(/3;i?) < p + (/3;R). For /i € 1, we define p((3,p) as the smallest density 
p of Gibbs measures P € "Pg G(/3, p); equivalently, as the left derivative, with 
respect to p, of the pressure p((3, p) = sup p (pp — /(/?, p)). Set 

p m := liminf p(/3, -m). 

Theorem 3.6. Let i? > n, m € (0,M) and i? ln > Vrf + 3 f m > Vd + 3r . 

(1) Suppose that v satisfies Assumption^ Then 

-v* < liminf/?^ 1 log p-{P;R). 

In addition, for every fixed v > v* , sufficiently large j3, p = exp(— fiis), 
there is a unique measure P in Qg(/3, p). It has no infinite cluster, P -almost 
surely, and satisfies 

P°(\C u (0)\ = k) = k Pk (P) < Cpexp(-f3ck) (15) 

for suitable C, c > and all k gN. 

(2) Suppose that v is continuous in (ro,ri).Then 

lim sup p + (p;Rm) < p m - 

13— »oo 

The constants C and c in Eq. (| 15[) can be chosen uniform in regions of the form 

P>p e ,p<exp{-P{v* + e))- 

If in addition v satisfies Assumptions[4]and[5l then p m is larger than the preferred 
ground state density po> an d in particular, bounded away from zero; this follows 
from Theorem 3.2 in |J12j . Moreover, we expect that for every R > r l5 as j3 — s- oo, 

lim P^ 1 log p-((3;R) = -v*, lim p+(/3;R) = p , 

[3 — >oo p— ^oo 

and for p-((3;R) < p < p+(j3; R) and very large j3, there should be non-ergodic 
Gibbs measures with percolation probability strictly between and 1 . This is what 
happens for lattice gases (see Appendix |A|) . For continuum systems, we have no 
proof of this conjecture; we have, however, a result pointing in the right direction: 
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Proposition 3.7. Suppose that v satisfies Assumptions ^ and [5l Let R > n. 

There are /3o,Po,C > such that for all ft > /3q, all p < po and all P G Gg((3,p), 
the following holds: if p = cxp(—[3is) > cxp(—[3is*) 7 then 

*■ — ' v — v 

k=l 

Thus at densities above cxp(— fiis*), the fraction of particles in finite-size clusters 
is small. 

4. Topology on V and continuity properties 

In this section we specify the topology on V that we use, recall some continuity 
properties of the functionals to be studied, and explain why the variational problems 
considered in this article have minimizers. We follow GZ93, G94 . 

Let M.(SX) be the set of finite measures on (fi,.F). The topology T£ of local 
convergence on .A/f(f2) is defined as follows. Let C be the class of measurable 
functions / : Q — > R that are local and tame, i.e., / G C if and only if there is a 
Borel subset B C R d and a constant c > such that / is a function of ujb alone 
and for all uj G SI, |/(w)| < c(l + Nb{u)). Then T£ is the smallest topology with 
respect to which all maps of the form P M> J„ /(w)P(du;), / G £, are continuous. 

The following holds [GZMIIUM] : 

• P and Pg are closed subsets of A4(£l). We endow them with the traces of 
the topology T£. 

• The Palm measure Pg — > A4(£l), P n- P° is continuous. 

• The particle density Pg — >• K, P M> /o(P) is continuous. 

• The entropy density Pg — >• M U {— oo}, P M> S'(P) is affine, upper semi- 
continuous, and has superlevel sets {S > — c} that are compact and se- 
quentially compact. 

• The energy density Pg — > R U {oo}, P i-> P(P) is lower semi-continuous. 
Furthermore, the cluster densities are continuous: 

Lemma 4.1. Por every k G N, £/ie map Pg — > R, P t— > pk(P), is continuous. 

Proof. Let gk(oo) := 1(|C W (0)| = fc). The function is local and bounded, thus in 
particular, tame. Therefore, by definition of T£, P h- > J^g^dP is continuous. Since 
P ^ Pk (P) is the composition of the latter map with the continous map P n- P° , 
it follows that P i— > Pk(P) is continuous. □ 

Now we can easily check that the variational problem in Theorem 13.11 admits a 
minimizcr. 

Lemma 4.2. Fix (3, p > and (pfc)fegN G [0,oo) N . Let .4 C Pg &e £/ie set o/ 
measures satisfying p(P) = p and Pk(P) = Pfe, /or every k G N. If A ^ the 
function A 3 P H- U(P) — f3~ 1 S(P) has a minimizer. 

Proof. If t/(P) — f3~ 1 S(P) = oo for every P G .A, there is nothing to show. If 
£/(P) — (3~ 1 S(P) < oo for some P G .A, let (P n ) be a minimizing sequence. The 
sequence (U{P n ) — /3 _1 S'(P„))„ e N is bounded from above and, because U(P) is 
bounded from below, S(P n ) is bounded from below too. Since the superlevel sets 
{S > — c} are r^-scqucntially compact, there is a subsequence P n converging to 
some P G Pg. The continuity of the maps p(-) and Pfc(-) ensures that P £ A, and 
the lower semi- continuity of U and —S shows that P is a minimizer. □ 
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5. Proof of Theorem 13.11 

The main idea for the proof of Theorem 13.11 is to apply a large deviations prin- 
ciple for the stationary empirical field proven in [G94, GZ93] and the contraction 
principle |DZ981 Section 4.2.1]. Two complications stand in our way. First, the 
large deviations principle in |G94( IGZ93] was shown in the grand-canonical rather 
than the canonical ensemble. Second, the cluster size densities can only be ex- 
pressed as functions of the stationary empirical field if we modify their definition 
and, loosely speaking, define them with periodic boundary conditions; this yields 
a modified variable p A CI . In order to circumvent these difficulties, we proceed as 
follows: 

• We show first that the large deviations principle in the canonical ensemble 
for p A implies a large deviations principle in the grand-canonical ensemble 
(Lemma 15. 1| . 

• We apply the contraction principle and show that /o A cr satisfies a large 
deviations principle with convex rate function (Lemma 15. 2| . 

• Next we show that (truncations of) p A and p A ° r , in the grand-canonical 
ensemble, are exponentially equivalent |DZ98[ Section 4.2.2]; this follows 
from Ruche's supcrstability bounds |R70j . As a consequence, the grand- 
canonical rate functions for p A and p A er must be equal (Lemma 15. 3[) . 

• Taking Legendre transforms, we deduce the desired identity for the canon- 
ical rate function (Lemma 15. 4|) . 

For the purpose of this section it is most convenient to work with measures that 
are not normalized and to suppress the /3-dependence in the notation. Let 

Qn%(A) := i J exp(~/3U{x 1 ,...,x N ))dx 1 ---dx N 

be a measure on A N with total mass Z\([3, N). The image of Q c ^\ under the map 
A N — > fl, . . . ,xn) h-> {xi, . . . ,xjv} is a measure on (f2, J 7 ), for which we use 
the same letter Q^A- Furthermore define 

oo 

Q^.a = 5 + J2 zN Qn% z = <*pG8aO- 

JV=1 

We use the same letter for the measure on Q and the measure on disjoint unions 
Un>qA n . A is a dummy space corresponding to u> — {0}: the event that there 
is no point at all has Q^A-measure 1. Remember that rf5 is equipped with the 
product topology and the corresponding Borcl tr-algebra. 

Lemma 5.1. Under Q^.a o,s |A| — > oo, N — > oo at fixed j3 > and the 
cluster size distribution p A := (pfe. A )fceN satisfies a large deviations principle with 
speed /3|A| and rate function 

We recall that f((3, •, •) : R + x M 1 ^ — > R is a lower semi-continuous, convex 
function, defined in all of p > and not only p S (0,p cp ). When p > p cp or 
Sfe+i kpk > P, it takes the value oo [JKMllj . 

Proof of Lemma \5.1l We consider Eq. (|16p as the definition of a function Jp^ and 
show that is a rate function for (pk,A)- Set z := exp(f3(j,). 
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Lower bound. Let O C R+ be an open set For every N and A, 

|tl A((Pfc 1 A)k eO)> z N Q%%({p k )k e O). 

Fix p > and apply the previous inequality to TV, A with A~ — > oo, |A| — > oo such 
that AT/|A| p. We get 

liminf TJTaT l °gQ»A(Pk,A)k S 0) > pp- inf /(/3,p, (p fe ))- 

P| A | (Pfc)6C 

Since the inequality holds for every p > 0, we can take the supremum over p > 
on the right-hand side, and to conclude note 

supfpp- inf /(p\p, (pfc))) =sup sup (HP- f(/3,p,(pk))) 

p>0 \ (Pfc)6C / p>0( Pfi ) e e> V 7 



sup sup 

(Pfc)ee>p>° 



(pp-/(p\p, (Pk)j) 



Upper bound. We use ideas explained in [R69| Section 3.4.5]. Let .4 C be a 
closed set. For every p > as |A| — > oo and AT/|A| —> p, 

limsup -L lng^X((p fc , A ) fc G A) < - inf /(p>, (p*))- (17) 
p|A| (p*)e>l 

Let b < oo be the stability constant from Eq. [TJ When Af/|A| — > 0, note that 

Q C N\((Pk,A) £A)< Z A ((3,N) < cxp(/3iV&)^! 
from which we get 

N 



limsup — logQ^ an A ((p fe ) £A)< limsup — ipNb-Nlog 
|A|-k» l A l |a|-k» |A| V |A|e 

Thus Eq. p7|) holds for all p > 0, provided we read the right-hand side as when 
p = 0. Fix po > 0. We claim that Eq. (|T7|) holds uniformly in p g [0, po]. More 
precisely, for every e > 0, there is a 5 > such that: for all p G [0, po] an d all AT, A 
satisfying |A| > 1/5 and |AT/|A| — p\ < 5, we have 

L ■ log QTA((pkA)k & A) < e— inf p, (p fe )). (18) 



p|A| (Pk)eA 

Indeed, if this was not the case, we could find e > 0, sequences (Nj), (Aj) and a 
p G [0,p ] such that |Aj| > j, |A^/|Aj | - p| < 1/j and 

lnQ C AfV((Pfe,A) fc e-4) >e - , inf /(p\p,(p fc )), 
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contradicting Eq. (|17[) : this proves the claim. Now because of the uniformity, for 
every e > and sufficiently large |A|, 

1 Lpo|A|J 
— log( E z N Q NA (( Pk , A ) G A) 
"I I iV=l 

< log(2 + PolA|) + e + , N _ n 

p|A| at V |A| (p h )e^l / 

< o(l) + e + sup( pp - inf /(/3, p, (p fc )) ) 

p>o^ (pfc)e-4 / 

= o(l) + e- inf Jp^Hpk))- 

If (pfc) = ^ A, the <$0 term corresponding, formally, to N = 0, docs not contribute 
to Q fj,,\{{pk,K) € .A). If G A the Jg-tcrm contributes a term 1 inside the loga- 
rithm, and we want to check that lim|A| _1 logl < — inf ^4 J/3. p. To this aim note 
that 

inf JpAiPk)) < mf (f(J3, P, 0) - ftp) < liminf (f(/3, p, 0) - pp) = 0. 

■4 p>0 p— W 

The last equality is shown by choosing a connected reference configuration {x\, . . . , x^-) 
which is such that (i) for suitable r > 0, every x G LlTves N x k=i ^( x w(k)^ r ) * s con " 
nected, (ii) on this set the energy is upper bounded by CN for suitable C > 0, 
and (iii) the balls are disjoint. Integrating in the neighborhood of x° given by 
xfL 1 B(xi, r) (and permutations of this set) yields the bound 



-Pf{fi, p,0)>-/3plog(|S(0,r)|e 1 - 



The right-hand side goes to as p — > 0. Thus we need not worry about the 
contribution from N = 0. 

It remains to estimate the terms from A" > po|A|. Choose po large enough so 
that zexp(^6+ l)//9 < 1/2. Remember AM > (N/e) N . Then 

E *"GS&(faM0 e A) < J2 z N ^™ 

N>p \A\ JV>p |A| 

/z|A|exp(/36 + 1)\ N , , 

^ E T "J <2cxp(-p |A|ln2). 



JV>Po|A| 



If inf^i < 00, we can choose /?o large enough so that po log 2 > /3 inf .4 J/3,p, and 
we obtain 

limsup-— rlogQ^A^A) G A) < - inf Jp !fl ((p k ))- (19) 
p|A| (Pfc)e^l 

If inf_4 = 00, we have /(/3,p, (pk)) = 00 for all p > and (p^) G A. But this 
implies that Q c ^ n A ((pk,A)k E A) = for all AT, A. Indeed, suppose by contradiction 
that there is a JV e N and a box A = [0, L] d such that Q c ^ n A ((pkA)k) e A) > 0. For 
neN, let L„ := n(L + R), A„ := [0,i„] d and A^ := n d N. Then 

«A„((Pfc,A) G A)) > (Q% n A ({pkA) G -4)) 
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and 

-p M.J(/3,p,( Pk )) > - {L l R)d logQ c ^ A ((Pfc,A) eA)> -oo, 

contradiction. Thus Q c ^ n A ((pk.A)k & A) = for all AT and A, whence Qfj,,A((pk.\) € 
.A) = 0; Eq. (IS) holds trivially. □ 

Next we define the stationary empirical field and the modified cluster size dis- 
tribution p A er . For uj e fl and A = [0, let w A cl := U fcez d6>Lfc(tJ n A) be the 
periodic continuation of cj n A. The translation invariant empirical field is 

ra,u = j-q- ^ i^p-dx. 

The Palm measure of i?A,w is 

Remark. The stationary empirical field associates with every configuration w a 
probability measure supported on configurations with the same relative coordinates 
as ui H A, but randomized center of mass. For example, in dimension d = 1, if u> 
consists of two particles 0, e, then R\ iUJ describes configurations {x, x + e} with x 
having uniform Lebesgue density 1/L. 

We define /c-clustcr densities with periodic boundary conditions as 
Mb"M : = ^l(|Ca(0)| = fc)P A)W (do>) = ^|{* e A | \C„ r (x)\ = k} 
We have = Pk(RA,w); compare with Eq. 0. Set p A (w) := ((p^ A (w)) feeN . 

Lemma 5.2. Under Q^a, P A ° r satisfies a large deviations principle with rate func- 
tion 

J Z(M = ^iU(P) - p-'SiP) - pp(P) | P e Ve, Vfe e N : Pfc (P) = 

Proof. Under Q^a, the stationary empirical field satisfies a large deviations prin- 
ciple with rate function U(P) — (3~ 1 S(P) — pp(P) |G94j . with compact sublevel 
sets. Since (p^.° A (uj))keN can be written as a continuous function of the stationary 
empirical field, we can apply the contraction principle |DZ98[ Section 4.2.1] and the 
result follows. □ 

Lemma 5.3. For every [3 > and p £ M, we have Jp^ = Jg „. 

Proof. When x £ A has distance > kR to the boundary <9A, we have C UA (a;) = 
C w p«(x). Set <9 fcjR A := {.t G A | dist(x,<9A) < fci?}. Then 

li. / \ r per, \ I , |wna fefl A| 

|«Pfc,AH - kp^aHI < — pq — 

and for every e > 0, 

p M)A (|p fc ,AH-p^H| > c ) < p„, a (aw > NA|). 

Here P^a : = (£a(/3, m)) ~ 1 Qh,a is the grand-canonical probability measure. Be- 
cause of the superstability of v, we know |R70j that for some suitable £ > such 
that all n-point correlation functions p\(x\, . ■ ■ , x„) are bounded in absolute value 
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by £ n ; £ depends on f3 and p but is independent of A or n. Integrating, we get 
bounds for factorial moments as 



E 



fi, A 



No kBA (N dkRA -!)••■ (N 9kRA - m + 1) < (£|0 fcfl A|) 



Let G(z) = J2 n =o P»,A-(Na kR (u)A = n)z n be the probability generating function of 
Nd kR A{u)- A Taylor expansion around z = 1 yields 



fc=0 



z=0 

A- 



n! ^ (ft - n)P j " n! ^ (k - n)\ 

k=n v ; k=n v y 

(£|d fe «A|)" fcflA| < ( e$\d kR k\ y „ dkRA 



For M g N such that M > 2e£|<9 fc ijA|, we have 

P, iA (^ BAM >M)<2(^^) M e^^. 
Applying the inequality to A/ = fce|A| we obtain, for sufficiently large |A| = L d 1 

P„ A (kAM - > e) < 2 (^) ktL \^-\ 

It follows that for every fixed K g N and e > 0, 

1 K 
lAKoo |Ap logP ^.A(Xl I/>*,a(w) - P P Ca(u)\ > ej = -oo, 

and the same identity holds with P^a replaced with Q^a (note Ha(/3, /i) > 1). 
As a consequence, the random variables (pi,A, ■ ■ • , Pk,a) and (pj a> ■ • ■ > Pk\) are 
exponentially equivalent in the sense of [DZ981 Section 4.2.1]. The contraction 
principle together with Lemma 15.21 and results from [JKMllj tell us that they 
satisfy large deviations principles; because of exponential equivalence, their rate 
functions must be equal |DZ98i Theorem 4.2.3]. The Dawson-Gartner theorem 
|DZ98[ Theorem 4.6.1] allows us to recover the rate function of the full vectors 
(k g N) from the rate functions of the truncated vectors (1 < k < K) and it follows 
that Jf}^ = J|™. □ 

Lemma 5.4. For every ft > 0, p S (0, p cp ) and (pk) € K+, we have 

f(/3,p, (p k )) = mi{U(P) - p- 1 S{P) \PeVe, p{P) = p, VfceN: p k (P) = p k }. 
Proof of Lemma\5^4\ Fix /3 > and {pk)keN & R+. For p > 0, set 

h( P ) := M{U(P) - pr l S{P) \PeVg, P (P) = p, Vfc g N : Pk (P) = Pk } 

and g(p) := f(f3, p, (p k )) for p g (0, p cp ), g(p) := oo for p < p cp ); the value of g(p cp ) 
can be determined in such a way that g is lower semi-continuous and convex in 
(0,oo). Lemma IS~3l implies that h and g have the same Legendre transforms, 

V/x g E : sup(p/x — = sup(/0/i — g(p)). 

p>0 p>0 

Since the functions h and g are lower semi-continuous in (0, oo) and convex [G94| 
IJKMllj . it follows that h(p) = g{p) for all p g (0, p cp ). □ 
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Proof of Theorem \3.1\ Thcorcm l3.1l is an immediate consequene of Lemmas 15 .41 and 

ma □ 



6. Proof of Theorem 13.21 and Prop. 

Proof of Theorem \3.2i (2)=> (1). Suppose that there is a shift-invariant Gibbs mea- 
sure P with p{P) — p and Pk{P) = Pk, for all k. By Theorem 13. 1[ 

min/(ftp,-) < /(ft p, ( Pk j) < U(P)-0- 1 S(P). 

On the other hand, let (p' k ) with Y^kLi kp'k — P- By Theorem 13.11 we find P' G Ve 
such that p(P') = p, p k (P') = p' k for all k, and f(/3,p, (p' k )) = U(P') - /T^P'). 
The Gibbs variational principle implies 

m Pi ( P ' k )) -w = u(p') - r'sip') - pp(p') > u(p) - p-'stp) - pp(p). 

Minimization over (p' k ) yields (remember p{P) = p) 

mmf(/3,p,-)>U(P)-p- 1 S(P). 

It follows that min/(ftp,-) = /(ftp, (p/0) = U{P) - ^"^(P), and (p k ) is a 
minimizer of /(ft p, •). 

(1)^(2). Conversely, let (p k ) be a minimizer of /(ftp, •). Thus /(ftp, (p*:)) = 
f(/3,p). By Theorem 13.11 we can find P G Vg such that p(P) = p, Pk{P) = Pk 
for all k, and /(ftp) = /(ftp, (p fc )) = C/(P) - ^" X 5(P). Wc know that /(ftp) = 
sup M (^p — p{f3,p)), and for every p G (0,p cp ), the supremum is actually attained. 
(This is because /i 4 £>(ft /«) is concave and continuous in R with limits 

— oo as p — > ±oo; note lim^oo d fl p(f3, p) = p cp .) Thus let p G K such that 
/(ftp) = p/x -p(ft/x). We get C/(P) - /^(P) - pp(P) = -p(^p) and as a 
consequence of the Gibbs variational principle, we find that P G Q((3,p). □ 

Proof of Prop, \3.3l (2) => f<?) By definition of our probability space f2, for every 
configuration u> and every bounded set A the number of particles in A is finite. 
Therefore, if u; has a cluster with infinitely many particles, this cluster has infinite 
diameter. 

(3) => fi?) If in a configuration to there is a cluster with infinite diameter, then this 
cluster must contain infinitely many particles - otherwise it would have diameter 
bounded by R times the cardinality of the cluster. 

(1 ) ^ (2) Let 

n 

g(w) := 1(|<U0)| = oo) = 1 - Urn V 1(|C W (0) = k). 

n— >oo * — ' 

fe=l 

The function is the limit of local, measurable functions and therefore measurable. 
Eq. (JSJ) gives 

I g(u>)P°(du>) = f l(|C(aOI = oo)P(dw). 

n n xeun[o,i] d 

Thus 

OO „ 

p{P) -Y^k Pk {P) = {xG^n[0,l] d | \\C u {x)\ = oo} P(du) (20) 
k=i Jn 

is the expected number of particles in [0, l] d belonging to an infinite cluster. If 
P(P) — J2kLi kpk{P) > 0, it follows right away that with positive probability, there 
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is an infinite cluster. If p(P) — Y)^-i kpk(P) = 0, using shift-invariance, we see 
that for every unit cube C(k), fc € 1 d (see p. [3]) the probability that the cube 
intercepts an infinite cluster is zero; it follows that the probability that there is an 
infinite cluster vanishes. □ 

7. Percolation properties of Gibbs measures 

In this section we prove Theorems 13.51 and 13.61 and Proposition 13.71 The proofs 
are a combination of results from |M751 |Z08, PY09, JKMTT] . what is known from 
cluster expansions IR69| IR70j , and equivalence of ensembles as in |G95j . 

7.1. Grand-canonical ensemble. For the proof of percolation, it is convenient to 
discretize space. Fix£>0. For k e Z d , let C(k) := [kit, (ki+l)£) x • ■ ■ x [k N £, (k N + 
1)£). Let C be the collection of cubes C(k), k e Z d . We say that two cubes are 
nearest neighbors if their centers have Euclidean distance i. A collection 1Z C C 
of cubes is connected if any two cubes in 7Z can be joined by a path (C\, . . . , C„) 
of cubes in 1Z such that for every j, the cubes Cj and Cj+i are nearest neighbors. 
The following lemma is a variant of well-known contour arguments |R69| Section 
5.3]. 

Lemma 7.1. Let P be a probability measure on {0, 1} . There is a constant ad > 
such that the following holds: if for some a > a e i, some A > 0, and all connected 
subsets 1Z C C . 

P(VC E 1Z : oj c = 0) < Aexp(-a|ft|), 
then P-almost surely, there is an infinite connected set W C C such that ujc = 1 
for all C e W. 

We refer to cubes C with ojc = as empty, and a cube with ujc = 1 as occupied. 

Proof of Lemma \7J\ Let T be the collection of the (d— l)-dimcnsional closed faces 
of the cubes in C. For example, [O,^]^" 1 x {0} e J. A set T C J is a contour if 
M. d \ UFerF splits into exactly two connected components, one finite (inside) and 
one infinite (outside). Let IntT C C be the collection of cubes that are inside T, 
and 9intr C IntT the cubes in the interior of T that touch the contour (sharing a 
face, an edge or a corner), i.e., 

a mt r = {C7eintr|<?n(U p)^0}. 

Note that <9i n tF is connected and |9mtF| > c|r|, for a suitable T-independent con- 
stant c. 

Fix Co € C. The number of contours T such that Co S IntT and |T| = n can 
be bounded by c\n d ^ d ~ 1 ' > exp(c2n), for some d-dependent constants c\,C2 > 0, see 
R69, Section 5.3]. It follows that if a is sufficiently large, 

oo 

P(lF : \T\ =n, Co G Intr, |<9 int r| is empty) < oo. 

71=0 

The Borcl-Cantelli lemma shows that P-almost surely, Co is enclosed in only finitely 
many contours with empty boundary 9j n tF. Thus we can pick a cell adjacent from 
outside to the union of all such contours. This cell cannot be empty, and it cannot 
be surrounded by another contour with empty boundary. It follows that it must be 
a member of an infinite connected set of occupied cells. □ 
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Proof of (2) in Theorem \3.5[ Pechersky and Yambartsev [PY09| proved a similar 
statement in dimension d = 2. We show that their proof can be adapted to d > 2. 
Let (ro, ri) be as in Assumption[TJ — m > inf/ r0)ri N v =: —M and f m g (ro, ri) such 
that v(f m ) < —to. Since v is continuous in (ro, ri), we can find ^, e, 5 > such that 
f m = I - e, < 2e < I, and v(r) < — m + 8 for all r g - 2e, £ + 2e]. Because of 
?m > ro, we also know that u(r) < for all r > £ — e. 

We partition R d into the cubes C(fc) of side-length I, as described before Lemma 
17.11 Let 7Z C C be a finite connected collection of cubes. Set n := \1Z\. Consider 
the events Q°(1Z) that there is no particle in := Uc^nC and ^(IZ) that every 
cube C £ K contains exactly one particle, and this particle has distance < e to the 
center of the cube. 

Choose L g £N large enough so that C [-L, L] d =: A and let Pa = Pp,n,&\$ be 
the grand-canonical (/?, /i)-Gibbs measure in A with free boundary conditions; write 
Sa = Sa|0(/3,/x) for the associated partition function. Fix an arbitrary numbering 
C\ , . . . , C n of the cells of 1Z and write u m for the center of C m . Note that if 1Z is 
connected, each cube has at least one neighbor. We have 



-i °° 7 n+k r r 
A £Z£ Kl JB(u\e) JB( 



dXr, 



u' L ,e ) 



x / dyi • • • dy fc e-^ xl -- x - l ' 1 '-' , '*\ 

./(A\A TC )* 

Because of our choice of £ and e, omitting the interaction between y^s and Xi's in 
the integrand can only decrease the Boltzmann weight. Since in addition for every 
x g x?P«e) 

J7(xi, ...,Xn) < (-m + S)(n - 1), 

we obtain 

Pa^CR.)) > z n \B(0,e)\ n e^ m - s ^ n -^P A (n Q (U)). 

The same argument applies to finite volume Gibbs measures with tempered bound- 
ary conditions £. Eq. ([3]) then shows that the previous inequality holds for every 
P g G(/3, n). It follows that, for every P g Q(fi, fi), 

P(Q (K)) < exp[-/3|ft|(/i + (1 - -L)(m - 8) - /T 1 log \B(0, t 

Lemma 17.11 applies and shows that if 

H > -m + S-/3- 1 log \B(0, e)| + fi~ x a d (21) 

then, P-almost surely, there is an infinite connected set of occupied cubes. Since the 
maximum distance between particles in two adjacent cubes is yd + 3 1, it follows 
that the particle configuration uj has P-almost surely an infinite P-cluster, for every 
R > Vd + 3£ = \/d + 3(r m + (-)■ 

To conclude, fix R > \J d + 3 f m . Choose 8, e, I as above, with the additional 
requirement that R > y/d~+3£. Let /3(<5, e) > (ad — log \B(0, e)\)/8. Then, for every 
> P(8, e), every fi > —m + 28, and every P g Q(P, p), there is an infinite cluster, 
P-almost surely. It follows that for sufficiently large f3, fi+(f3;R m ) < —m + 28. 
Since 8 can be chosen arbitrarily small, the proof is complete. □ 
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Proof of (1) in Theorem \3.5[ For f3 > and k G N, define the cluster partition 
functions 

Zf{[3) := ^ [ e-P u( > - X2 <-- Xk h({0,X2T..,x k } is connected) dx 2 ... dx k 



><nfc-l 



and let i? cl (/3) be the radius of convergence of J^SL, z ^kiP): compare |M75[ Prop. 
3.3]). Let P G C?(/3,/x) be a Gibbs measure that can be obtained as a limit of 
finite volume, grand canonical Gibbs measures with empty boundary conditions. 
Suppose that z < R cl (f3). Then P(there is an infinite cluster) = [M75| Theorem 
3.1]. Miirmann's proof moreover yields, for every cube A C M. d and all fceff, the 
bound 



P[3x 6wnA: \C u (x)\ = k^ 

<fr / e -^ xi --' x "h(3j G : Xj €A)dx 1 ...dx k 

k\ J( R <i)fc 

< k\A\z k Zl\f3). 

Let £ > be as a in the proof of Lemma 15.31 If P is shift-invariant, Eq. ([5]) and 
Rucllc's superstability bounds |R70] show that 

|A|P°(|C u (0)| = k) = f {x€unA:\C u (x)\=k}\)P(duj) 

= P(3x G u n A : \C u {x)\ = fc) + 0((£|A|) 2 ) 

as |A| — > 0. It follows that 

Pk {P) < z k Zl\f3). (22) 
In order to go from shift-invariant limits of finite volume Gibbs measures to general 
Gibbs measures, we use the theory of Mayer expansions. It is well-known that 
for every f3 > 0, some strictly positive i? May (/3) > 0, and activities z = exp(/3 / u) < 
R May (f3), there is a unique Gibbs measure Pp,n G G(f3, At); furthermore, the pressure 
and the correlation functions admit absolutely convergent expansions in powers of 
z (with temperature-dependent coefficients) |R701 Theorem 5.7]. The measure Pp^ 
is shift-invariant. Since every finite volume Gibbs measure converges to an infinite 
volume Gibbs measure, this limit must be Pp^, and we can apply the previous 
considerations to Pp^- Thus we see that when z < iiim(i? cl (/3), i? May (/?)), there 
is a unique Gibbs measure. It is shift-invariant, has no infinite cluster (P-almost 
surely), and satisfies Eq. ([22]) . 
Next, recall 

Zf{(i) < e-P ke "e k \B(0,R)\ k -\ i? May (/3) > ^rrrjp 

where |||w||| := |P(0,r hc )| + f< x < >rh \v(\x\)\dx, see [JKM111 Proof of Lemma 4.1] or 
[M75l Proof of Prop. 3.1] for the first inequality and jPU09| Theorem 2.1] for the 
second inequality. We deduce 

liminf/r 1 log min (i? cl (/?), i? May (^)) > eoo . 

Since every /i_ < fi^ 1 logmin(i? cl (/3), R May (f3)) satisfies condition (fTTj) . we obtain 
the desired lower bound on /i-(/3; R). The upper bound (fT2|) of Pk(P) follows from 
Eq. d22j) and the upper bound of Zf(f3). □ 
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7.2. Canonical ensemble. Theorem 13.61 is deduced from Theorem 13.51 with the 
help of a good control of the density as a function of the chemical potential. Recall 
that the pressure p(/3, p) = sup p (p/i — /(/3, p)) is a convex function of p. Therefore 
the derivative p(/3,p) := cLp(/3, p) exists almost everywhere and is an increasing 
function of p. Moreover, if p = p(f3,p), then Qg((3 7 p) C Q(f3 7 p). If p(f3,p) has 
different left and right derivatives p\ and p T with respect to /x, the previous inclusion 
holds for every p £ [p\, p r ] . 

Proof of (1) in Theorem \S.b\ Remember i? May (/3) and i? cl (/3) from the proof of (1) 
in Theorem 13.51 For cxp(/3^) < i? May (/3), the pressure is differentiablc in p and 
p(f3,p) is well-defined. Fix e > and let f3 e > large enough so that for f3 > f3 e , 
cxp(^[ eoo - e]) < min(i? May (/3),i? cl (/3)). For [3 > j3 e and p < p(f3, eoc - e), we 
know that p = p(f3,p) for a unique p; for this p, Qe{P,p) = Q{P,p) = {P} with a 
unique Gibbs measure P. From the proof of Theorem I3.5( we see that P assigns 
probability zero to the event that there is an infinite cluster. As a consequence, 
P-(f3; R) > p{(3, eoo - e). Since 

lim /T 1 log p{(3, eoo - e) = - inf (e h - k(e x - e)) (23) 

fl^oo k£N\ / 

and inf e>0 mik^n[Ek — k(e 00 — e)] = v* |J12j . we deduce p~(f3; R) > exp(— f3v*(l + 
o(l)))- ' 

In addition, if v > v* and (3 is sufficiently large, we see that p = p(/3, p) for 
a unique p < /3" 1 log min(i? May (/3), R cl (/?)), and the remaining part of Theorem 
13.61 follows from the corresponding statement in Theorem 13.51 The claim on the 
uniformity of the constants stated after Theorem 13.61 follows by combining the 
inequalities |JT5J|, ||22J), and fl23J|. □ 

Proof of (2) in Theorem Fix e > 0. By Theorem [3~5l there is a f3 e > such 
that for every (3 > fl e , p+{(3; R m ) < —m + e. Let p(/3, —m + e) be the left derivative 
of p{f},p) with respect to p at p = —m/2 + e. If p > p e (/3), then Qg(f3,p) C 
Q({3,p) for some p > —m + e. When f3 > f3 e , it follows that this p is larger 
than p + ((3;R m ), and Theorem 13.51 tells us that for every P € Ge{(3,p) C Q(f3,p), 
there is an infinite cluster, P-almost surely. Thus p + ((3;R m ) < p(f3,—m + e) and 
liminf/^oo p+(/3, /u) < p m . □ 

Proof of Prop. \3. 7\ Prop. 13.71 is a consequence of Theorem 1.8 in [JKMllj and 
our Theorem 13.21 the result from jJKMllj says that the bound of our proposition 
holds for every minimizer (pk)ken of /(/?, p, ■), and Theorem 13.21 says that for every 
P € Qe{(3, p), the vector (pk(P))keN is a minimizer of /(/3, p, •). □ 

Appendix A. Lattice gas and Ising model 

Here we explain how our results relate to known results for lattice systems. 
We shall be very sketchy and refer the reader to to the review [GHMOlj and the 
references therein for precise definitions and proofs. 

Consider the nearest neighbor lattice gas in Z d (d > 2) with hard-core on- 
site interaction and attractive nearest-neighbor interaction. As is well-known, this 
model can be recast as an Ising model. Occupied lattice sites (n x = 1) arc mapped 
to spin a x — +1 and empty lattice sites (n x = 0) to spin a x = — 1; thus a x = 2n x — l. 
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Formally, we have 

J ^ — > ^ — > J ^ — > dJ ~\~ p ^ — > 
-— 2^ n xn v - ^2_^n x = ~— °xO y ^ — > y + const, 

\x-y\ = \ x \x-y\ = \ x 

which determines the external magnetic field h of the Ising model in terms of the 
chemical potential p of the lattice gas as h = (p + dJ)/2. The parameter J > 
measures both the strength of attraction between particles and the strength of the 
ferromagnetic coupling between spins. 

The relevant notion of percolation is dependent site percolation. Dependent refers 
to the underlying measure on {0, 1} Z , which is a Gibbs measure rather than a prod- 
uct of Bernoulli measures, and site percolation refers to the notion of connectivity 
- two occupied lattice sites x, y £ II 1 are connected if there is a path (x±, . . . , x n ) of 
nearest neighbor (|xj+i — Xj\ = 1) joining x\ = x and x n = y, such that each site of 
the path is occupied, n Xj = 1. In the Ising picture, we are interested in percolation 
of + clusters. 

Our first remark is that at low temperature, the phase transition and the per- 
colation transition coincide; this observation, as alluded to in the introduction, is 
the driving motivation for the present article's investigation. Fix a temperature 
T > and vary the external field h (or the chemical potential p) . At temperatures 
above the Curie temperature Tq > 0, there is no phase transition (the pressure 
stays analytic and the Gibbs measure is unique); at T < Tq, there is a first-order 
phase transition at h = 0. On the other hand, let P^ h be the Gibbs measure with 
+ boundary conditions. We know that for all T > 0, there is a threshold h(T) E E 
such that the P^h -probability of having an infinite cluster is at external fields 
h < h(T), and 1 at fields h > h(T). Furthermore, there is a temperature T + > 
such that h(T) = for T <T + and h{T) < for T > T+ [ABL87] . and it is known 
that T + < Tq in high dimensions. Thus at temperatures above the Curie temper- 
ature, there is a percolation transition but no phase transition; at temperatures 
between T + and Tq, as h is increased, the percolation transition happens before the 
phase transition; and below T + , the percolation transition coincides with the phase 
transition. 

The next observation is that the value at which the transition takes place is 
consistent with our grand-canonical Theorem 13.51 indeed, the zero external field 
h = corresponds, in the lattice gas picture, to a chemical potential \x = —dJ, 
which can be interpreted as a ground state energy per particle - if in Z d every 
lattice site is occupied, the energy per particle is eoo = — dJ. 

Finally, this agreement of thresholds extends to the canonical ensemble. At T < 
Tq and external field h = 0, there are two shift-invariant Gibbs measures h , with 
magnetizations m±(T). Low temperature contour expansions show that as T — > 0, 
m±(T) = ±l+0(exp(-J/T)). Transforming as p±(T) = (m±(T) + l)/2, we obtain 
that the corresponding curves for the lattice gas satisfy p-(T) = 0(exp(— J/T)) and 
p + (T) = 1 + 0(cxp(- J/T)), which compares nicely with TheoremES for T < T+, 
the magnetizations m±(T) of the Ising model are also percolation thresholds (this 
is not true for T + < T < T c ). 
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